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In this paper we present a summary of results obtained for scalar field theories 
using the Feynman-Schwinger (FSR) approach. Specifically, scalar QED and x 2( P 
theories are considered. The motivation behind the applications discussed in this 
paper is to use the FSR method as a rigorous tool for testing the quality of com- 
monly used approximations in field theory. Exact calculations in a quenched theory 
are presented for one-, two-, and three-body bound states. Results obtained indi- 
cate that some of the commonly used approximations, such as Bethe-Salpeter ladder 
summation for bound states and the rainbow summation for one body problems, pro- 
duce significantly different results from those obtained from the FSR approach. We 
find that more accurate results can be obtained using other, simpler, approximation 
schemes. 



I. INTRODUCTION 

In the study of hadronic physics one has 
to face the problem of determining the quan- 
tum dynamical properties of physical sys- 
tems in which the interaction between the 
constituents is of a nonperturbative nature. 
In particular, such systems support bound 
states and clearly nonperturbative methods 
are needed to describe their properties. As- 
suming that such systems can be described by 
a field theory, one has to rely on some approx- 
imation scheme. One common approxima- 
tion is known as perturbation theory. Pertur- 
bation theory involves making an expansion 
in the coupling strength of the interaction. 
The Green's function in field theory can be 
expanded in powers of the coupling strength. 
In order to be able to obtain a bound state 
result one must sum the interactions to all or- 
ders. Most practical calculations to date have 
been done within the Bethe-Salpeter frame- 
work, where the resulting kernel is perturba- 



tively truncated. In this paper we will dis- 
cuss another method, which is based on the 
path integral formulation of Feynman and 
Schwinger [l[ H| • It was initiated by Simonov 
and collaborators @, H H H, 0, II in their 
study of quantum chromo-dynamics (QCD). 

With the discovery of QCD, nonperturba- 
tive calculations in field theory have become 
even more essential. It is known that the 
building blocks of matter, quarks and glu- 
ons, only exist in bound states. Therefore 
any reaction that involves quarks will nec- 
essarily involve bound states in the initial 
and/or final states. This implies that even 
at high momentum transfers, where QCD 
is perturbative, formation of quarks into a 
bound state necessitates a nonperturbative 
treatment. Therefore it is essential to develop 
new methods for doing nonperturbative cal- 
culations in field theory. 

The plan of this article is as follows. In the 
following section a brief review of the parti- 
cle trajectory method in field theory will be 
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given. The Feynman-Schwinger representa- 
tion will be introduced through applications 
to scalar fields. In particular, the emphasis 
will be on comparing various nonperturba- 
tive results obtained by different methods. 
It will be shown with examples that non- 
perturbative calculations are interesting and 
exact nonperturbative results could signifi- 
cantly differ from those obtained by approx- 
imate nonperturbative methods. 

Results for scalar quantum electrodynam- 
ics are discussed in section 3. In particular, 
a cancellation is found between the vertex 
corrections and self-energy contributions, so 
that the exact result can be described by es- 
sentially the sum of only the generalized ex- 
change ladders. In section 4 the bound states 
are obtained for a x 2 4> theory in a quenched 
approximation. Exact results are presented 
for the binding energies of two- and three- 
body bound states and compared with rela- 
tivistic quasi-potential predictions. Section 5 
deals with the stability of the x 2( P theory. It 
is argued that the quenched theory does not 
suffer from the well known instability of a <f) 3 
theory. The paper closes with some conclud- 
ing remarks. 



II. FEYNMAN-SCHWINGER 
REPRESENTATION 

Nonperturbative calculations can be di- 
vided into two general categories: (i) inte- 
gral equations, and (ii) path integrals. In- 



tegral equations have been used for a long 
time to sum interactions to all orders with 
various approximations fiol 11 . H2I k|. In 
general a complete solution of field theory to 
all orders can be provided by an infinite set of 
integral equations relating vertices and prop- 
agators of the theory to each other. How- 
ever solving an infinite set of equations is be- 
yond our reach and usually integral equations 
are truncated by various assumptions about 
the interaction kernels and vertices. The 
most commonly used integral equations are 
those that deal with few-body problems. The 
Bethe-Salpeter equations [lfj] are the start- 
ing point of those investigations. Approxima- 
tions schemes have extensively been studied, 
where in addition to solutions of the Bethe- 
Salpeter equation [13, ll^ , f\ i 
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dimensional reductions [1 
have been explored for the N-particle free 
particle Greens' function. In most calcula- 
tions the ladder approximation has been used 
for the kernel of the resulting equations. Is- 
sues of convergence of these schemes remain. 
Therefore it is important to have exact solu- 
tions of field theory models available to test 
these approximations. One promising way to 
reconstruct exact solutions of field theory is 
the path integral method. 



Path integrals provide a systematic 
method for summing interactions to all or- 
ders. The Green's function in field theory is 
given by the path integral expression: 
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While path integrals provide a compact ex- gral is a nontrivial task. 

pression for the exact nonperturbative result In general, field theoretical path inte- 

for propagators, evaluation of the path inte- grals must be evaluated by numerical inte- 
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gration methods, such as Monte-Carlo in- 
tegration. The best known numerical in- 
tegration method is lattice gauge theory. 
Lattice gauge theory involves a discretiza- 
tion of space-time and numerical integrations 
over field configurations are carried out us- 
ing Monte-Carlo techniques. A more effi- 
cient method of performing path integrals in 
field theory has been proposed and it con- 
sists of explicitly integrating out the fields. 
It has been demonstrated to be highly suc- 
cessful for the case of simple scalar interac- 
tions 0, M, H H S H H M, 0. 

This method is known as Feynman-Schwinger 
representation (FSR). Through applications 
of the FSR, the importance of exact nonper- 
turbative calculations will be shown with ex- 
plicit examples. 

The basic idea behind the FSR approach 
is to transform the field theoretical path inte- 
gral (0) into a quantum mechanical path in- 
tegral over particle trajectories. When writ- 
ten in terms of trajectories, the exact results 
decompose into separate parts, and permit 
us to study the individual role and numeri- 
cal size of exchange, self-energy, and vertex 
corrections. This, in turn, alows us to study 
different approximations to field theory, and, 
in some cases, prove new results. 

To illustrate these ideas, we consider the 
application of the FSR technique to scalar 
QED. The Minkowski metric expression for 
the scalar QED Lagrangian in Stueckelberg 
form is given by 



C 



sqed = -m\ 2 - ^F 2 + X -ix 2 A 2 

+(d^-ieA^)x*(d^+ieAnx 



-^(d-A) 2 , (2) 

where is the gauge field of mass /i, x is the 
charged field of mass m, = d^A" — d v A^ 
is the gauge field tensor, and, for example, 
A 2 = A^N 1 . The presence of a mass term 
for the exchange field breaks gauge invari- 
ance, and was introduced in order to avoid 
infrared singularities that arise when the the- 



ory is applied in 0+1 dimensions. For dimen- 
sions larger than n = 2 the infrared singu- 
larity does not exist and therefore the limit 
\i — > can be safely taken to insure gauge 
invariance. 

The path integral is to be performed in 
Euclidean metric. Therefore we perform a 
Wick rotation: 



exp 



d xC 
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exp 
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(3) 



The Wick rotation for coordinates is obtained 
by 



x — > -IXq, 



d 



d 

dx° 



id • 



(4) 



The transformation of field A under Wick ro- 
tation is found by noting that under a gauge 
transformation: 



A^ -> A^ + d^A. 



Then, under a Wick rotation, 
A -> iA , 



(5) 



(6) 



and the Wick rotated Lagrangian for SQED 
becomes 
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SQED — X 
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(7) 



The exchange field part of the Lagrangian is 
given by 

£. 4 = i/l^V - \d„d„)A p + if 2 



2 M 



fJ,U 



+(1 - \)dpd v 



A„ 



(8) 



We employ the Feynman gauge A = 1 which 
yields 



C A = -A u (fx 2 - D)A U . (9) 
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where 



Se — J d 4 X C-SQEDi 



(11) 



FIG. 1: The dashed lines represent exchanges of 
the gauge field with mass /x and the solid lines 
the propagation of the matter fields with mass 
m. Note the matter loop in one of the middle 
exchanges. All loops of this kind are neglected in 
the quenched approximation (when detS" =1). 

The two-body Green's function for the tran- 
sition from an initial state x) to final 
state Qf(y,y) is given by 



and a gauge invariant 2-body state $ is de- 
fined by 

$(x,x) = x*(x)U(x,x) X (x). (12) 

The gauge link U(x, y) which insures gauge 
invariance of bilinear product of fields is de- 
fined by 



U(x, y) = exp 



— ie / dzA(z) 



(13) 



One can easily see that under a local gauge 
transformation 



X{x) -► e—^'x(x) 



G(y,y\x,x) = N / V x * J V x J VA 

x^*^ e - SE , (10) 



J,eA(x), 

A^x) -> A^x) + d^A(x), (14) 



$i(x, x) remains gauge invariant 



$(x, x) — > exp 



-ieA{x) + ieA(x) — ie dz^ d^A 



i 



X*(x)U(x, x)x(x) = $(x,x). (15) 







Performing path integrals over x an d X* fields in Eq. (jlOjl one finds 

G(y,y\x,x) = N(vA(detS) U(x, x)U*(y, y)S(x, y)S(x, y) e 



-S[A] 



(16) 



where the interacting 1-body propagator 
S(x, y) is defined by 



S(x,y) = (y 



m? + H(z,p) 



with 



H(z,p) = (p + ieA(z)) 2 . 



The Green's function Eq. (116)1 includes con- 
tributions coming from all possible interac- 
tions. The determinant in Eq. (|Trj|) accounts 



for all matter (xx) loops. Setting this de- 
terminant equal to unity (detS* — > 1, referred 
to as the quenched approximation) eliminates 
x) (17) all contributions from these loops (illustrated 
in Fig. and greatly simplifies the calcula- 
tion. 

Analytical calculation of the path integral 
over the gauge field A in Eq. (fTrTj) seems dif- 
ficult due to the nontrivial A dependence in 
S(x,y). In more complicated theories, such 
as QCD, integration of the gauge field, as 
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Z(X) y 



FIG. 2: The contour C, known Wilson loop, 
that arises in Eq. (|24|). 

far as we know, cannot be done analytically. 
Therefore, in QCD, the only option is to do 
the gauge field path integral by using a brute 
force method. Here one usually introduces 
a discrete space-time lattice, and integrates 
over the values of the field components at 
each lattice site . However, for the simple 
scalar QED interaction under consideration, 
it is in fact possible to go further and elim- 
inate the path integral over the field A. In 



order to be able to carry out the remaining 
path integral over the exchange field A it is 
desirable to represent the interacting prop- 
agator in the form of an exponential. This 
can be achieved by using a Feynman repre- 
sentation for the interacting propagator. The 
first step involves the exponentiation of the 
denominator in Eq. (|17|): 

POO 

S(x,y)= / dse- sm2 (y\exp[-sH}\x) . (19) 
Jo 

This expression is similar to a quantum me- 
chanical propagator with s = it and H a 
Hamiltonian which is a covariant function 
of 4-vector momenta and coordinates. It is 
known how to represent a quantum mechani- 
cal propagator as a path integral. The repre- 
sentation is in terms of the Lagrangian, and a 
covariant Lagrangian can easily be obtained 
from the Hamiltonian ()18j) 



H(z,p) = (p + ieA(z)f 



i 2 



L(z, z) — — — iez ■ A(z) . 



(20) 



Using this Lagrangian, the path integral representation for the interacting propagator be- 
comes 



S{x,y) 



(Vz) xy exp 



-sm 



dr s z (t s ) -ie dr s zA(z(r s )) 



(21) 



where Zi(r s ) is a particle trajectory which is a parametric function of the parameter r s , with 
s > t s > and endpoints Zi(0) = Xi, Zi(s) = yi, and i — 1 to 4. This representation 
allows one to perform the remaining path integral over the exchange field A. The final 
result for the two-body propagator involves a quantum mechanical path integral that sums 
up contributions coming from all possible trajectories of the two charged particles 



G 



ds 



| ds | (Vz) xy | (Vz) 



xy 



-K[z,s]-K[z,s] 



(W(C)), 



where the parameter r s is rescaled, so that t s = sr, the kinetic term K is defined by 

1 f 1 

K[z,s] = m 2 s-\ / drz 2 (r), 

and the Wilson loop average (W(C)) is given by 

(W(C)) = J VAexp -ie j dz A(z) -~J d^zA{z){ii 2 - d 2 )A(z 



(22) 



(23) 



(24) 
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where the contour of integration C (shown in Fig. |2J) follows a clockwise trajectory x — > 
y^y^x-^xas parameters r, and f are varied from to 1. The A integration in 
the definition of the Wilson loop average is of standard gaussian form and can be easily 
performed to obtain 



(W(C)) 



cxp 



e' 
~~2 



dz„ 



c 



dz v A^y (z - Z,/l) 



c 



D ipx 



(27r)> 2 + /i 2 



When it is necessary to regulate the ultraviolet singularities in 
subtraction will be used, so that (J26|) will be replaced by 



d*p e^(A?-/i 2 )(A^-/i 2 ) 



(25) 
(26) 

a double Pauli-Villars 



A^(x,/i) = g^ 



}(27r) 4 (p 2 + /i 2 )(p 2 + A 2 )(p 2 + A|) 



(27) 



Through the results given in Eqs. (J22j) . (J23), 
and either (|2l)Jl or (J2ZJ), the path integration 
expression involving fields has been trans- 
formed into a path integral representation in- 
volving trajectories of particles. 

Equation (|22j) has a very nice physical in- 
terpretation. The term A^ u (z a — Zb, fi) de- 
scribes the propagation of gauge field iter- 
ations between any two points on the parti- 
cle trajectories, and the appearance of these 
interaction terms in the exponent means 
that the interactions are summed to all or- 
ders with arbitrary ordering of the points 
on the trajectories. Self-interactions come 
from terms with the two points z a and Zb on 
the same trajectory, generalized ladder ex- 
changes arise if the two points are on differ- 
ent trajectories, and vertex corrections arise 
from a combination of the two. Because the 
particles forming the two-body bound state 
carry opposite charges, it follows that the self 
energy and exchange contributions have dif- 
ferent signs. 

The bound state spectrum can be deter- 
mined from the spectral decomposition of the 
two-body Green's function 



G{T) = J2 



-m n T 



(2* 



n=0 



tween the initial and final states 



T = -(y 4 + y i -x i - x 4 ) 



(29) 



where T is defined as the average time be- 



In the limit of large T, the ground state mass 
is given by 

r d JVZS'[Z]e- s W 

(30) 



III. SCALAR QUANTUM 
ELECTRODYNAMICS 

In this section we take a closer look at 1- 
body mass pole calculations for the case of 
SQED. Two popular methods frequently used 
to find the dressed mass of a particle are to do 
a simple bubble summation, or to solve the 
1-body Dyson-Schwinger equation in rainbow 
approximation. It is interesting to compare 
results given by the bubble summation and 
the Dyson-Schwinger with the exact FSR re- 
sult. Below we first give a quick overview 
of how dressed masses can be obtained in 
bubble summation and the Dyson-Schwinger 
equation approaches. For technical simplic- 
ity, the 1-body discussion will be limited to 
0+1 dimension. 
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The simple bubble summation involves a 
summation of all bubble diagrams to all or- 
ders. The dressed propagator is given by 



Bubble sum: 



A,(p) 



p2 _|_ m 2 _|_ 



(31) 



The dressed mass M is determined from the 
self energy using 



M = ^Jm 2 + Z(iM). 



(32) 



The self energy for the simple bubble sum (in 
0+ld) is given by 



E(p) 



1 



dk 



(2p - kf 
[{p — k) 2 + m 2 ] 



1 . (33) 



The self energy integral in this case is trivial 
and can be performed analytically, and the 
dressed mass is determined from Eq. ()32)) 

The rainbow Dyson-Schwinger equation 
sums more diagrams than the simple bubble 
summation (Fig. |3J). The self energy of the 
rainbow Dyson-Schwinger equation involves 
a momentum dependent mass: 



E(p) 



-! 



dk 



1 



x 



_ 0O 2n(k 2 + ^) 

(2p - A;) 2 
- A;) 2 + m 2 + E(fc)] 



- 1 . (34) 



In this case the self energy is nontrivial and it 
must be determined by a numerical solution 
of Eq. (J34|) . The dressed mass is determined 
by the logarithmic derivative of the dressed 
propagator in coordinate space 

M = ^log[A d (t)]. (35) 

The type of diagrams summed by each 
method is shown in Fig. El Note that the 
matter loops do not give any contribution 
as explained earlier. Results obtained by 
these three methods are shown in Fig. HJ 
It is interesting to note that the simple 



\ i 
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Dyson-Schwinger: 



♦ • •• 



Feynman-Schwinger: 



4 • • • • • — 4-^ — 4 — k- 



FIG. 3: Various interactions included in each 
approach are shown. 



FSR 

- Bubble sum 
Dyson-Schwinger 
Bubble (unphys.) 




FIG. 4: The function M(e 2 ), calculated in 0+ld 
with values of m = fi = 1 GeV, using the FSR 
approach, the Dyson-Schwinger equation in the 
rainbow approximation, and the bubble summa- 
tion. While the exact result is always real, the 
rainbow DSE and the bubble summation results 
become complex beyond a critical coupling. 

bubble summation and the rainbow Dyson- 
Schwinger results display similar behavior. 
While the exact result provided by the FSR 
linearly increases for all coupling strengths, 
both the simple bubble summation and the 
rainbow Dyson-Schwinger results come to a 
critical point beyond which solutions for the 
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r(p.q) 




FIG. 5: Because of the interactions, the one- 
particle irreducable vertex functions T n (n = 
3, 4, • • • ) depend on the external momenta. 

dressed masses become complex. This ex- 
ample very clearly shows that conclusions 
about the mass poles of propagators based 
on approximate methods such as the rainbow 
Dyson- Schwinger equation can be misleading. 

In general a consistent treatment of 
any nonperturbative calculation must involve 
summation of all possible vertex corrections. 
Vertex corrections are those irreducible dia- 
grams that surround an interaction vertex. 
The elementary vertex is the three-point ver- 
tex, T 3 , but the particle interactions will lead 
to the appearance of nth order irreducible 
vertices, r„, as illustrated in Fig. The 
propagation of a bound state therefore in- 
volves a summation of all diagrams with the 
inclusion of higher order vertices (Fig. |BJ). A 
rigorous determination of all of these vertices 
is not feasible. In the literature on bound 
states r n>3 interaction vertices are usually 
completely ignored. The 3-point vertex T 3 
can be approximately calculated in the lad- 
der approximation |40||. However a rigorous 
determination of the exact form of the 3-point 
vertex is not possible, for this requires the 
knowledge of even higher order vertices. 

In order to be able to make a connec- 
tion between the exact theory and predictions 
based on approximate bound state equa- 
tions it is essential that the role of interac- 
tion vertices be understood. The Feynman- 
Schwinger Representation (FSR) is a useful 
technique for this purpose. 




r 

5 



FIG. 6: Exact computation of the two-body 
bound state propagator requires the summation 
of all particle self energies, vertex corrections, 
and ladder and crossed ladder exchanges. 

We now present the interesting outcome, 
that the full bound state result dictated 
by a Lagrangian can be obtained by sum- 
ming only generalized ladder diagrams ( "gen- 
eralized" ladders include both crossed lad- 
ders and, in theories with an elementary 
four-point interaction, both overlapping and 
non-overlapping "triangle" and "bubble" di- 
agrams) . 

We adopt the following procedure for de- 
termining the contribution of vertex correc- 
tions in 3+1 dimension. We start with an 
initial bare mass m and calculate the full 
two-body bound state result with the in- 
clusion of all interactions: generalized lad- 
ders, self energies and vertex corrections. Let 
us denote the result for the exact two-body 
bound state mass by M2 0t (e 2 ,m), since it 
will be a function of the coupling strength 
e and the bare input mass m, and the su- 
perscript "tot" implies that all interactions 
are summed. Next we calculate the dressed 
one-body mass Mi(e 2 ,m). Then using the 
dressed mass value Mi(e 2 ,m) we calculate 
the bound state mass Mf xch (e 2 ,Mi) by sum- 
ming only the generalized exchange interac- 
tion contributions. In this last calculation we 
sum only exchange interactions (generalized 
ladders) , but the self energy is approximately 
taken into account since we use the (con- 
stant) dressed one-body mass as input. How- 
ever the vertex corrections and wavefunction 
renormalization are completely left out since 
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C(p,q) 




p-q 



S(p,m(p)) 



S(pfm(p)) 



e 2 /(47t) 

FIG. 7: Two-body bound state mass for SQED 
in 3+1 dimensions. Solid triangles are 2Mi, 
open squares are Mf xch , and open triangles are 
M^ ot . Here p/m = 0.15, Ai//i = 3, and A 2 /fi = 
5. The smooth lines are fits to the "data" . Note 
that M| xch = Mf ot to within errors. 



we use the original vertex provided by the 
Lagrangian. In order to compare the full re- 
sult where all interactions have been summed 
with the result obtained by two dressed par- 
ticles interacting only by generalized ladder 
exchanges we plot the bound state masses ob- 
tained by these methods. Numerical results 
are presented in Fig. [7| This result is quali- 
tatively similar to that obtained analytically 
for SQED in 0+1 dimension (33[. 

The numerical results presented here yield 
the following prescription for bound state cal- 
culations: In order to get the full result for 
bound states it is a good approximation to 
first solve for dressed one-body masses ex- 
actly (summing all generalized rainbow dia- 
grams), and then use these dressed masses 
and the bare interaction vertex provided by 
the Lagrangian to calculate the bound state 
mass by summing only generalized ladder in- 
teractions (leaving out vertex corrections). In 
terms of Feynman graphs this prescription 
can be expressed as in Fig. |S] 

The significance of the results presented 
above rests in the fact that the problem 
of calculating exact results for bound state 



r:(p,q.k) 



p-q-k 



S(p,M) I S o (p-q,M) 
1 



p £ p-q-k 
— a • 

/ \ 
q'' \ 



/ \ 

q ' "k 



FIG. 8: The correct two-body result can be 
obtained by simply using a dressed constituent 
mass and a bare vertex, and ignoring the contri- 
butions of higher order vertices. 

masses in SQED has been reduced to that of 
calculating only generalized ladders. Sum- 
mation of generalized ladders can be ad- 
dressed within the context of bound state 
equations [13, EH- Here we have shown 
the connection between the full prediction 
of a Lagrangian and the summation of gen- 
eralized ladder diagrams. Our results are 
rigouous for SQED, but are only suggestive 
for more general theories with spin or internal 
symmetries. Since we have neglected charged 
particle loops (our results are in quenched ap- 
proximation) , and the current is conserved in 
SQED, it is perhaps not surprising that the 
bare coupling is not renormalized, but the 
fact that the momentum dependence of the 
dressed mass and vertex corrections seem to 
cancel is unexpected. If we were to unquench 
our calculation, or to use a theory without a 
conserved current, it is reasonable to expect 
that both the bare interaction and the mass 
would be renormalized. 

Finally, we call attention to a remarkable 
cancellation that occurs in the one-body cal- 
culations. The exact self energy shown in 
Fig. [7| (and also in Fig. 0] for different param- 
eters) is nearly linear in e 2 j2£| • This remark- 
able fact implies that the exact self energy is 
well approximated by the lowest order result 
from perturbation theory. It is instructive to 
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see how this comes about. If we expand the 
self energy to fourth order, expanding each 
term about the bare mass m, we have 

S d \p 2 ) = m 2 -p 2 + Y>(p 2 ) 

= m 2 - p 2 + S 2 + (p 2 - m 2 )S' 2 
+£ 4 , (36) 

where = S^(m 2 ) is the contribution of or- 
der e e evaluated at p 2 = m 2 , £' = d£(p 2 ) / dp 2 
evaluated at p 2 = m 2 , and the formula is 
valid for p 2 —m 2 ~ e 2 . Expanding the dressed 
mass in a power series in e 2 

M 2 = m 2 + m 2 2 + ml^ , (37) 

where m 2 is the contribution of order e , and 
substituting into Eq. (|3b|h give 

Ml = m 2 + S 2 + S 2 S 2 + E 4 + • • • ,(38) 

The mass is then 

2m 8m 3 
+ ■••, (39) 



The linearity of the exact result implies that 
the forth order term in Eq. ()39|) must be zero 
(or very small), and this can be easily con- 
firmed by direct calculations! 

The cancellation of the fourth order mass 
correction (and all higher orders) is remi- 
nisent of the cancellations between general- 
ized ladders that explains why quasipotential 
equations are more effective that the ladder 
Bethe-Salpeter equation in explaining 2-body 
binding energies. It shows that a simple eval- 
uation of the second order self energy at the 
bare mass point is more accurate than solu- 
tion of the Dyson Schwinger equation in rain- 
bow approximation. 

The general lesson seems to be that at- 
tempts to sum a small subclass of diagrams 
exactly is often less accurate than the ap- 
proximate summation of a larger class of di- 
agrams. 

In the next section we consider the appli- 
cation of the FSR approach to scalar x 2 in- 
teraction 



IV. SCALAR x 2 <P INTERACTION WITH THE FSR APPROACH 

We consider the theory of charged scalar particles x °f mass m interacting through the 
exchange of a neutral scalar particle of mass //. The Euclidean Lagrangian for this theory 
is given by 



(40) 



The 2-body propagator for the transition from the initial state = X*( x )x( x ) to final state 
= X*(y)x(y) is given by 



G(y,y\x,x) = N 



^Vx*^Vx^V(f) $}$iexp - J d 4 xC 



(41) 



After the usual integration of matter fields is done, the Green's function reduces to 



G(y,y\x,x) = N Dtp (detS) S(x,y)S(x,y)e 



-L [ 



(42) 
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with the free Lagrangian, L , and the interacting propagator, S(x,y), defined by 



L [<f>] = - I #z <Hz){f - %)4>(z) 



S{x,y) = (y 



1 



m? + H(z,p) 



x) 



with the Hamiltonian 



H(z,p) = p 2 - g4>(z) 



(43) 



(44) 



As in the case of scalar QED we employ the quenched approximation: det5* — > 1. 

We exponentiate the denominator by introducing an s integration along the imaginary 
axis with an e prescription 



S(x,y) 



rioo 

/ ds e - s{m2+ie) (y | exp[-sH] \ x) 
Jo 



(45) 



This representation should be compared with the representation used earlier in SQED 
Eq. ([19)1 . Here the integration is done along the imaginary axis because H is not posi- 
tive definite. 

Again, a quantum mechanical path integral representation can be constructed by recog- 
nizing that Lagrangian corresponding to the H of Eq. (|44|) is given by 



L{z,z) = — + g<f){z) . 
The path integral representation for the interacting propagator is therefore 



S(x,y) 



—i 



T>z exp 



4 



is{m 2 + it) / dTz 2 (r)+ig / dr 4>{z{t)) 



G = -J ds J ds J (Vz) xy | {Vz) m e ^,s]+iK[z,s] Iif 



where mass and kinetic term is given by 



1 r 1 

K[z,s] = {m 2 + ie)s / drz 2 (r). 

4s Jo 



The field integration is a standard gaussian integration 



dT<P(z(r))+ / df<f>{z{f))) -L 



= J 1)0 exp +ig(^ 
= exp(-V [z, s)-2 V 12 [z, z, s, s) - V [z, s]j . 



(46) 



(47) 



This representation allows the elimination of the integral over the exchange field <fi. The 
2-body propagator reduces to 



(48) 



(49) 



(50) 
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where V and Vi 2 (self and exchange energy contributions in Fig. EJ) are defined by 



V [z,s] = ^— s 2 dr\ dr'A(z(r) -z(r'),fi) 



s v\z\ 



Vi 2 [z,z,s,s] 



2 /l /l 

— ss dr 
2 Jo Jo 



dr A(z(t) — z(r),fj,) . 



(51) 



(52) 



It should be noted that the interaction terms explicitly depend on the s variable, which was 
not the case for SQED. The interaction kernel A is given by 



A(x,/i) 



I 



d 4 p 



(27r)V + /i 2 



An 2 \x 



■Ki(n\x\) . 



(53) 



In order to be able to compute the path integral over trajectories, a discretization of the 
path integral is needed 



(Vz) 



XII 



(N/4ns 



,2N 



d Zi . 



(54) 



The s dependence is crucial for correct normalization. After discretization, the 1-body 
propagator takes the following form 



G 



2N 



ds 



im 2 s — % s 2 v\z] 

4s 



(55) 



where v[z] was defined in Eq. f!51|) . This is an oscillatory and regular integral and it is 
not convenient for Monte-Carlo integration. The origin of the oscillation is the fact that s 
integral was defined along the imaginary axis, 



Rep. 1 : S(x,y) =< y 



ds exp 



—s(m 2 — d 2 + gcj) + ie) 



x > 



In earlier works 0, H3] a nonoscillatory Feynman-Schwinger representation was used, 



Rep. 2 : S(x,y) =< y\ dsexp 



I. 



-s(m 2 - d 2 + gcj)) 



x > 



Rep. 2 leads to a nonoscillatory and divergent result 

° ds 



G oc 



-2N 



exp 



k 2 

2 2 r n 

m s h s v\z\ 

4s 



(56) 



(57) 



(58) 



and the large s divergence was regulated by a cut-off A. This is not a satisfactory prescription 
since it relies on an arbitrary cut-off. Later it was shown [28|, |30| that the correct procedure 
is to start with Rep.l and make a Wick rotation such that the final result is nonoscillatory 
and regular. The implementation of Wick rotation however is nontrivial. Consider the 
s-dependent part of the integral for the 1-body propagator 



G oc 



1, 



ds 

l2lv ex P 



im s 



k 2 



s 2 v \z] 



(59) 



It is clear that a replacement of s — > is leads to a divergent result. The problem with the 
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Im(s) 




Re(s) 

FIG. 10: Wick rotation in the s integration. 

Wick rotation (Fig ITUj) comes from the fact 
that the s integral is infinite both along the 
imaginary axis and along the contour at in- 
finity. These two infinities cancel to yield a 
finite integral along the real axis. As g — > 
the dominant contribution to the s integral 
in Eq. ()59jl comes from the stationary point 



is 



k 

2m 



(60) 



Therefore one might suppress the integrand 
away from the stationary point by introduc- 
ing a damping factor R 

R(s,s ) = l-(s-is ) 2 /T 2 . (61) 

With this factor the integral in Eq. (foTJj) is 
modified 



Got 



OO 7 

as 

* 



im s 



.k 2 s 2 v[z] 
'As R 2 (s,s Q )_ 



(62) 



This modification allows us to make a Wick 
rotation since the contribution of the contour 



FIG. 11: As the effective coupling strength g 2 is 
increased, the stationary point disappears. 

at infinity now vanishes. However this pro- 
cedure relies on the fact that there exists a 
stationary point. It can be seen from the 
original expression Eq. (|59|) that this is not al- 
ways true. According to the original integral 
the stationary point is given by the following 
equation 



im + i 



Introducing s 



4s 2 



/s 



2m 



2sv[z] 
and g 2 



this equation becomes 
-1 + s 2 



gV 



o. 



. (63) 
k (v[z\) /m 3 

(64) 



The stationary point is determined by the 
first intersection of a cubic plot with the pos- 
itive s axis as shown in Fig. ^2 The plot in 
Fig. ITTI shows that, as the effective coupling 
strength g 2 is increased the curve no longer 
crosses the positive s axis. Therefore be- 
yond a critical coupling strength the station- 
ary point vanishes and mass results should be 
unstable. Limiting discussion to cases where 
the original expression Eq. ()59|) has a crit- 
ical point, we now turn to perform a Wick 
rotation on the modified expression Eq. [)62jl . 
The Wick rotation in Eq. (|62|) amounts to 
a simple replacement s — > is, and a regular, 
nonos dilatory integral is found: 



G(x 



OO 7 

as 

5 



n k 2 s 2 v\z] 
-m 2 s h 



4s R 2 (is,s ) 



(65) 
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FIG. 12: Number of steps a particle takes be- 
tween initial and final coordinates is discretized. 
The space-time is continuous and there are no 
space-time boundaries 

At first look it seems that the new integral 
always has a stationary point determined by 
the following equation 

9 k 2 2sv\z] 
-m 2 + — + 1 

4s 2 R 2 (is, s ) 



— s 2 v[z 



R 2 (is,s ) 



0. (66) 



The key point to remember is that the sta- 
tionary point we find after the Wick rotation 
should be the same stationary point we had 
before the Wick rotation. This is required to 
make sure that the physics remains the same 
after the Wick rotation. Self consistency 
therefore requires that the stationary point 
after the Wick rotation is at s = isq. In that 
case R(is , s ) = 1, and (R 2 (is , s ))' = 
and Eq. (JBTtj) determining the critical point 
reduces to the earlier original Eq. (JfiUj) . 

The regularization of the ultraviolet sin- 
gularities is done using Pauli-Villars regular- 
ization, which is particularly convenient for 
numerical integration since it only involves a 
change in the interaction kernel 

A(x,fi) — ► A(x, fi) — A(x, a/i). (67) 

Calculations of the x 2 4> interaction in 
3+ld require numerical Monte-Carlo integra- 
tion. The first step is to represent the parti- 
cle trajectories by a discrete number of iV + 1 
points with boundary conditions given by 



Z = X = (Xi, 22,^3,0) 

z N = y= (yi,y2,V3,T) ■ 



(6S 



going from the initial time to the final time 
along a trajectory in a 4-d coordinate space. 
This is very different from the discretization 
employed in lattice gauge theory. Contrary 
to lattice gauge theory, in the FSR approach 
space-time is continuous and rotational sym- 
metry is respected. An additional, impor- 
tant benefit is the lack of space-time lattice 
boundaries, which allows calculations of arbi- 
trarily large systems using the FSR approach. 
This feature provides an opportunity for do- 
ing complex applications such as calculation 
of the form factors of large systems. 

In doing Monte- Carlo sampling we sam- 
ple trajectories (lines) rather than gauge field 
configurations (in a volume). This leads to 
a significant reduction in the numerical cost. 
The ground state mass of the Green's func- 
tion is obtained using 



m 



fVZS'[Z)e 



-S[Z] 



,-S\Z] 



(69) 



Sampling of trajectories is done using the 
standard Metropolis algorithm, which insures 
that configurations sampled are distributed 
according to the weight e~ s ^ z \ In sampling 
trajectories the final state (spacial) coordi- 
nates of particles can be integrated over, 
which puts the system at rest and projects 
out the 5-wave ground state. As trajectories 
of particles are sampled, the wave function of 
the system can be determined simply by stor- 
ing the final state configurations of particles 
in a histogram. 

In sampling trajectories the first step is 
thermalization. In order to insure that the 
initial configuration of trajectories has no ef- 
fect on results, the first 1000 or so updates are 
not taken into account. Statistical indepen- 
dence of subsequent samplings is measured 
by the correlation function X(n), defined as 



X(n) 



(m(i)m(i + n)) — (mY 
(m) 2 



(70) 



The discretization employed in the FSR is for where m(i) is the mass measurement at the 
the number of time steps a particle takes in z'th update. 
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In order to insure that the location of 
the stationary point is self consistent, as dis- 
cussed earlier, its location must be deter- 
mined carefully. The stationary point can 
be parametrized by so = CT/(2m), where 
T/2m is the location of the stationary point 
when the coupling strength g goes to zero. 
As the coupling strength is increased, the 
stationary point moves to larger values of 
So (recall Fig. ITTj) and C increases. Even- 
tually a critical value of the coupling con- 
stant is reached beyond which there is no self 
consistent stationary point. In order to be 
able to do Monte-Carlo integrations, an ini- 
tial guess must be made for the location of 
the stationary point. Self consistency is re- 
alized by insuring that the peak location of 
the s distribution in the Monte- Carlo inte- 
gration agrees with the initial guess for the 
stationary point [30]. In Fig. [Till the depen- 
dence of the location of the stationary point 
on the coupling strength for 2-body bound 
states is shown. The figure shows that be- 
yond the critical point g 2 ~ 100 GeV 2 , C goes 
to infinity implying that there is no station- 
ary point. A similar critical behavior was also 
observed in Refs. (3?], EH within the context 
of a variational approach. In Fig. El FSR 2- 
body bound state mass results are shown for 
m x = 1 GeV, = 0.15 GeV. These results 
are all for a Pauli-Villars mass of 3/x. Also 
are shown the predictions of various integral 
equation calculations. 

The FSR calculation sums all ladder and 
crossed ladder diagrams, and excludes the 
self energy contributions. According to 
Fig. EH all bound state equations underbind. 
Among the manifestly covariant equations 
the Gross equation (labeled GR in Fig. H3J) 
gives the closest result to the exact calcula- 
tion obtained by the FSR method. This is 
due to the fact that in the limit of infinitely 
heavy-light systems the Gross equation effec- 
tively sums all ladder and crossed ladder dia- 
grams. The equal-time equation (labeled ET 
in Fig. IHJ) also produces a strong binding, 
but the inclusion of retardation effects pushes 
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FIG. 13: The dependence of the peak of s- 
distribution on the coupling strength for the 2- 
body bound state is shown. The peak location is 
given by sq = CT/2m. There is no real solution 
for C beyond g 2 = 100 GeV 2 . 
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FIG. 14: The coupling constant dependence of 
the 2-body bound state mass is shown. Beyond 
the critical coupling strength of g 2 = 100 GeV 2 
the 2-body mass becomes unstable. The Bethe- 
Salpeter equation in the ladder approximation 
gives the smallest binding. The other models 
are described in the text. 

the ET results away from the exact results 
(Mandelzweig- Wallace equation ji^l, labeled 
MW in Fig. [T3j). In particular the Bethe- 
Salpeter equation in the ladder approxima- 
tion (labeled BSE in Fig. IT4"j) gives the low- 
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FIG. 15: The dependence of the peak of s- 
distribution on the coupling strength for the 3- 

body bound state mass is shown. The peak lo- 

• • u r^^r/n °-ru • i FIG. 16: 3-body bound state results for 3 equal 

cation is given by so = C 1 /2m. 1 here is no real J H 

solution for C beyond g 2 = 81 GeV 2 . 



mass particles of mass 1 GeV. 



est binding. Similarly the Blankenbecler 



Sugar-Logunov-Tavkhelidze equation [20, 21 
(labeled BSLT in Fig. ITljl gives a very low 
binding. A comparison of the ladder Bethe- 
Salpeter, Gross, and the FSR results shows 
that the exchange of crossed ladder diagrams 
plays a crucial role. 

In Fig. HS1 the dependence of the loca- 
tion of the stationary point on the coupling 
strength for 3-body bound states is shown. 
In Fig. [TB] the 3-body bound state results for 
3 equal mass particles of mass 1 GeV are 
shown. For the 3-body case the only avail- 
able results are for the Schrodinger and Gross 
equations. According to the results presented 
in Fig. Unj the bound state equations under- 
bind for the 3-body case too. The Gross 
equation gives the closest result to the exact 
FSR result. 

Determination of the wavefunction of 
bound states is done by keeping the final 
state configurations of particles in a his- 
togram. For example, for a 3-body bound 
state system, the probability distribution of 
the third particle for a given configuration of 
first and second particles is shown in Fig. El 
In the upper-left panel of Fig. [T7[ the two 



fixed particles are very close to each other so 
that the third particle sees them as a single 
particle. However as the fixed particles are 
separated from each other the third particle 
starts having a nonzero probability of being 
in between the two fixed particles (separa- 
tion increases as we go from upper right to 
lower left panels in Fig. 117)1 . Eventually when 
the two fixed particles are far away from each 
other the third particle has a nonzero prob- 
ability distribution only at the origin (the 
lower left panel in Fig. 117]). 

Up to this point the FSR method has been 
derived and various applications to nonper- 
turbative problems have been presented. In 
the next section we discuss the stability of 
the x 2< P theory. 



V. 



STABILITY OF THE SCALAR X 2 <, 
INTERACTION 



Scalar field theories with a x^X'P interac- 
tion (which we will subsequently denote sim- 
ply by x 2( P) have been used frequently with- 
out any sign of instability, despite an argu- 
ment in 1952 by Dyson[35] suggesting insta- 
bility, and a proof in 1959 by G. Baym H~~ 
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FIG. 17: The four panels show the evolution of the probability distribution for the 3rd particle as 
the distance between the two fixed particles is increased. When the fixed particles are very close 
to each other the third particle sees them as a single particle (the upper left plot). As the fixed 
particles are separated, the third particle starts penetrating between them (2nd and 3rd panels), 
and when two fixed particles are far apart (as shown in the lower right panel), the third particle is 
most likely to be found between the two fixed particles. 



showing that the theory is unstable. For ex- 
ample, it is easy to show that, for a limited 
range of coupling values < g 2 < g 2 rit , the 
simple sum of bubble diagrams for the prop- 
agation of a single x particle leads to a stable 
ground state, and it is shown in Ref. |l3[ that 
a similar result also holds for the exact re- 
sult in "quenched" approximation. However, 
if the scalar x 2 4> interaction is unstable, then 
this instability should be observed even when 
the coupling strength g is vanishingly small 
g 2 — > + , as pointed out recently by Rosen- 
felder and Schreiber|37| (see also Ref. [39l|). 
Both the simple bubble summation and the 



quenched calculations do not exhibit this be- 
havior. Why do the simple bubble summa- 
tion and the exact quenched calculations pro- 
duce stable results for a finite range of cou- 
pling values? 

A clue to the answer is already provided 
by the simplest semiclasical estimate of the 
ground state energy. In this approximation 
the gound state energy is obtained by mini- 
mizing 

E = m 2 X 2 + ^ 2 <p 2 -g<Px 2 , (71) 

where m is the bare mass of the matter par- 
ticles, and /i the mass of the "exchanged" 
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quanta, which we will refer to as the mesons, expanded in terms of creation and annihila- 
The minimum occurs at tion operators that depend on time 



E = mY - g 2 ^ . (72) x ^ r) = J dk m [a{k) e' 1 ^ + b\k) e ik ' x ) 

This is identical to a x 4 theory with a cou- ^ r ) = f dk^ [c(k) e~ ik ' x + c\k) e ik ' x ] , 
pling of the wrong sign, as discussed by Zinn- J 

Justin The ground state is therefore sta- (74) 
ble (i.e. greater than zero) provided 

where x = {t, r} and 

9 <Scnt = — (73) _ d , k 

dkm " (2ir) 3 2E (k) ^ 
This simple estimate suggests that the the- y J ml ; 

ory is stable over a limited range of couplings 

if the strength of the x field is finite. We with E U k ) = Vm 2 + k 2 . The equal-time 

now develop this argument more precisely commutation relations are 
and show under what conditions it holds. 

Before presenting new results, we lay the [a{k), a^(k')] = (2tt) 3 2E m (k) 5 3 (k - k') . 

foundation using the variational principle. In (76) 

the Heisenberg representation the fields are The Lagrangian for the x 2< P theory is 



C = x t [ d 2 -m 2 + gcfi] X + \ <fi (d 2 - 4> , (77) 

and the Hamiltonian if is a normal ordered product of interacting (or dressed) fields (fi d and 
Xd 

H[4> d , Xd ,t] = I d A r : { ( + (y X df + m\l 



+ 2 



1 ^ 



dt 



+ (v<fi d y + fi 2 



-grAd}:. (78) 



This hamitonian conserves the difference between number of matter and the number of 
antimatter particles, which we denote by uq. Eigenstates of the Hamiltonian will therefore 
be denoted by |n , A), where A represents the other quantum numbers that define the state. 
Hence, allowing for the fact that the eigenvalue may depend on the time, 

H[<f> d)X d,t] \no,X) =M no , x (t) \no,X) • (79) 
In the absence of an exact solution of (|79|h we may estimate it from the equation 

M n0tX (t) = (n , X\H[(fi d ,Xd,t] |n , A) 

= (n , A| U~\t, 0)H[(f>, X , 0)U{t, 0) \n , A) 

= (no,X,t\H[(f),x,0]\no,X,t) , (80) 
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where U(t, 0) is the time translation operator which carries the Hamiltonian from time t = 
to later time t. We have also chosen t = to be the time at which the interaction is turned 
on, 4>d(t) = U^it, 0)(f)(0)U(t, 0), and the last step simplifies the discussion by permitting 
us to work with a Hamiltonian constructed from the free fields 4> and x- [If the interaction 
were turned on at some other time to, we would obtain the same result by absorbing the 
additional phases exp(±iEto) into the creation and annhilation operators.] 
At t — the Hamiltonian in normal order reduces to 



H[4>, X ,0] = J dk m E m (k)Af (k,k) + J dp^E^p)c ] {p)c(p) 



2 J uj(k-k') 

where 



c\k' -k) + c{k - k') 



(81) 



Af (k, k') = {a\k)a{k') + tf(k)b(k')} , 

J\fi(k, k') = AT {k, k') + {a\k)b\-k') + a(-k)b(k')} (82) 

and uj{k) = \J [i 2 + k 2 . To evaluate the matrix element ()80|) we express the the eigenstates 
as a sum of free particle states with n matter particles, n pair pairs of XX particles, and £ 
mesons: 



j oo oo 

\n 0l X,t) = \n ,a(t),(3(t)}} = -— ^ a„ pair (t)/3 £ (t) |n , n pair , I) , (83) 

™ ' ™ pair =0 £=0 

where 7(t) is a normalization constant (defined below), the time dependence of the states is 
contained in the time dependence of the coefficients a(t) and f3(t), and 

\n ,n pair ,£)= / === (84) 

•/ V( w + ™pair ) ^pair « *! 



with m = n + n pair , n 2 = n pair and 

/„ ni ?i2 ^ 

= / n ^/(^) n «*&•/(<&■) n (85) 

n 1 « 1 I 1 



i=l j = l 1=1 



The particle masses in dk and dp have been suppressed; their values should be clear from 
the context. The normalization of the functions f(p) and g(p) is chosen to be 

J dkf\k) = J dpg 2 (p) = l (86) 



which leads to the normalization 



(n , X,t | n , A,t) = 1 , (87) 
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if 7 (t) = a(t)p(t) with 

<x> 

a\t)= <^ (*)=«(*)• 

^pair 
OO 

/5 2 W= £ = /?(*)•/?(*) ■ (88) 

The expansion coefficients {ctn pair (£)} and {/^(i)} are vectors in infinite dimensional spaces. 

In principle the scalar cubic interaction in four dimensions requires ultraviolet regular- 
ization. However the issue of regularization and the question of stabilty are qualitatively 
unrelated. For example, the cubic interaction is also unstable in dimensions lower than four, 
where there is no need for regularization. The ultraviolet regularization would have an ef- 
fect on the behavior of functions f(p), and g(p), which are left unspecified in this discussion 
except for their normalization. 

The matrix element (|8U|) can now be evaluated. Assuming that f(k) = f(—k) and 
g{k) = g(—k), it becomes: 



M no> x{t) = {n + 2L(t)}m + G(t) jl - gV {n + 2L(t) + 2L x {t)} y/G^t) , (89) 
where the constants m, p, and V are 

rh = J dkE m (k) f(k) , p = J dp E^(p) g 2 (j>) 

v _ f dk m dk' m f(k)f(k')g(k-k') 
J v/m 2 + (k - k') 2 

and the time dependent quantities are 

L[t) - ^ — ' G{t) -l^y^) 



^pair 1 



a 2 (t) 



Note that L and G are the average number of matter pairs and mesons, respectively, in the 
intermediate state. 

The variational principle tells us that the correct mass must be equal to or larger than 
(JBHj). This inequality may be simplified by using the Schwarz inequality to place an upper 
limit on the quantities L\ and G\. Introducing the vectors 

fx = {ax,V2a 2 , • ■ ■} = {\/n a n } 

f 2 = {y/n + 1 a , y/n + 2 ax, ■ ■ ■} = Wn + n a n -i} 
h= {(3x,V2/3 2 ,---} = {Vi(3 e }, (92) 
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we may write 



Ll{t) ~ « 2 (t) - ^5) 

= ^L(t){n + l + L(t)} 



Hence, suppressing explicit reference to the time dependence of L and G, Eq. (|89j) can be 
written 



M n0iX {t) > (n + 2L)m + Gjl-gV ^n + l + L + VZy -1^ VG . 



(94) 



Minimization of the ground state energy with 
respect to the average number of mesons G 
occurs at 

VG~o = ^r {(v^o + l + ^ + v/T) 2 -lj . 

^ (95) 
At this minimum point the ground state en- 
ergy is bounded by 

M no . x {t)> {n + 2L}m-fiG . (96) 

If we continue with the minimization process 
we would obtain M no< \{t) — ► — oo as L — > oo, 
providing no lower bound and hence suggest- 
ing that the state is unstable. However, if 
L is finite , this result shows that the ground 
state is stable for couplings in the interval 
< g 2 < glu with 



2 

9crit 



4 jttm (no + 2L) 



2 ' 



(97) 

This interval is nonzero if the number of mat- 
ter particles, n , and the average number of 
XX pahs, L, is finite. In particular, if there 
are no Z diagrams or \X loops in the inter- 
mediate states, then the ground state will be 
stable for a limited range of values of the cou- 
pling. 

This result also suggests strongly that the 
system is unstable when g 2 > g 2 Tit , or when 



L — > oo (implying that g 2 Tit — > 0). However, 
since Eq. ()96|) is only a lower bound, our argu- 
ment does not provide a proof of these latter 
assertions. 

We now discuss the effect of the Z-graphs 
and of the matter loops on the stability. 
Using the Feynman-Schwinger representation 
(FSR), we will show that the ground state is 
(i) stable when Z-diagrams are included in 
intermediate states, but (ii) unstable when 
matter loops are included. 

The covariant trajectory z(t) of the par- 
ticle is parametrized in the FSR func- 
tion of the proper time r. In theory the 
FSR expression for the 1-body propagator for 
a dressed x-particle i n quenched approxima- 
tion in Euclidean space was given qualita- 
tively in Eq. (|65|). The detailed expression, 
needed in the following discussion, is 



G(x,y) = I ds 
'o 



N 
Atts 



2N N-l 



n * 



x exp^ -K[z, s] - V[z, s r ] \ , (98) 



where the integrations are over all possible 
particle trajectories (discretized into N seg- 
ments with N — 1 variables z$ and boundary 
conditions z = x, and zn = y) and the ki- 



22 



netic and self energy terms are 



N 



N 



K[z,s] = m 2 s + — ^(z i -^„ 1 ) 2 ,(99) 



i=l 



V[z,s] 



2 2 N 

9 s 2 



2N 2 



^A(^ ii)At ) , (100) 



where A(z,}i) is the Euclidean progagator 
of the meson (suitably regularized), Sz^ = 



Za ~t~ Zi—\ Zj Zj 



), and 



R(s,s ) i + (s - s ) 2 /r 2 ' 



(101) 



was 



(The need for the substitution s — > 
discussed above in Sec. IV.) 

In preparation for a discussion of the ef- 
fects of Z-diagrams and loops, we first dis- 
cuss the stability of Eq. (|9~H]) when neither 
Z-diagrams nor loops are present. To make 
the discussion explicit, consider the one body 
propagator in 0+1 dimension. Since the in- 
tegrals converge, we make the crude approxi- 
mation that each z% integral is approximated 
by one point (since we are excluding Z- 
diagrams, the points may lie along the clas- 
sical trajectory). If the boundary conditions 
are z = and z N = T the points along the 
classical trajectory are z { = iT/N, and 

AT N T 2 
T r r i 2 X ^ / \ 2 2 

K[z, s\=m s+— } \Zi—Zi-\) =ms+— . 



4s 
(102) 

If the interaction is zero, this has a stationary 
point at s = so = T/(2m), giving 



K[z, s] = Kq = mT . 



(103) 



yielding the expected free particle mass m. 
[Note that half of this result comes from the 
sum over (zj — Zi_i) 2 .} The potential term 
(|100j) may be similarily evaluated; it gives a 
negative contribution that reduces the mass. 

We now turn to a discussion of the effect 
of Z-diagrams. For the simple estimate of 
the kinetic energy, Eq. ()102j) . we chose inte- 
gration points Zi = iT/N uniformly spaced 
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x 

1 z(x) 
folded trajectory: 



• y 




FIG. 18: It is possible to create particle- 
antiparticle pairs using folded trajectories. How- 
ever folded trajectories are suppressed by the 
kinematics, as discussed in the text. 

along a line. The classical trajectory con- 
nects these points without doubling back, so 
that they increase monotonically with proper 
time, r. However, since the integration over 
each Zi is independent, there also exist tra- 
jectories where Zi does not increase mono- 
tonically with r. In fact, for every choice 
of integration points Zj there exist trajecto- 
ries with z-i monotonic in r and trajectories 
with Zi non-monotonic in r. The latter dou- 
ble back in time, and describe Z-diagrams in 
the path integral formalism. Two such tra- 
jectories that pass through the same points 
Zi are shown in Fig. These two trajecto- 
ries contain the same points, Zi, but ordered 
in different ways, and both occur in the path 
integral. 

Now, since the total self energy is the 
sum of potential contributions V[z, s] from 
all (z^ Zj) pairs, irrespective of how these co- 
ordinates are ordered, it must be the same 
for the straight trajectory z(t) and the folded 
trajectory Zf(r): 



V[Zf,8]=V[z,8]. 



(104) 

However, according to Eq. (|99|). the kinetic 
energy of the folded trajectory is larger than 
the kinetic energy of the straight trajectory 



K[z f ,s] > K[z,s] , 



(105) 



because it includes some terms with larger 
values of (zi — Zi-i) 2 . Since the kinetic en- 
ergy term is always positive, the folded tra- 
jectory (Z-graph) is always suppressed (has 



23 



continuous deformation: 




• y 



Z(T) 



FIG. 19: A folded trajectory at the end point of 
the path, and a similar one with Z\ closer to zq. 

a larger exponent) compared with a corre- 
sponding unfolded trajectory (provided, of 
course, that g 2 < g 2 ctit ). 

This argument holds only for cases where 
the trajectory does not double back to times 
before zq = or after zn = T '. An example 
of such a trajectory is shown in Fig. (up- 
per panel). Here we compare this folded tra- 
jectory to another folded trajectory, z'f, with 
point Z\ closer to the starting point z (lower 
panel of Fig. IT9|). This new folded trajectory 
has points spaced closer together, so that the 
kinetic energy is smaller and the potential en- 
ergy is larger, and therefore 

K[z f ,s]-V[z f ,s}>K[z' f ,s}-V[z' f ,s}. 

(106) 

It is clear that the larger the folding in the 
trajectory, the less energetically favorable is 
the path, and the most favorable path is 
again an unfolded trajectory with no points 
outside of the limits z < Z\ < z N . 

While these arguments have been stated 
in 0+1 dimensions for simplicity, they are not 
dependent on the number of dimensions, and 
can be extended to the realistic case of 1+3 
dimensions. 

We conclude that a calculation in 
quenched approximation, where the creation 
of particle-antiparticle pairs can only come 
from Z-graphs, must be more stable (pro- 
duce a larger mass) than a similar calcula- 
tion without any XX pairs. The quenched 
X 2 <f) theory therefore is bounded by the same 
limits given in Eq. ([97)1 . This conclusion sup- 
ports, and is supported by, the results of 



Refs. which show, in the quenched 

approximation, that the x 2( P interaction is 
stable for a finite range of coupling strengths. 

It is now clear that the instability of x 2 <P 
theory must be due to either (i) the possibil- 
ity of creating an infinite number of closed 
XX loops, or (ii) the presence of an infinite 
number of matter particles (as in an infinite 
medium). Indeed, the original proof given 
by Baym used the possibility of loop creation 
from the vacuum to prove that the vacuum 
was unstable. 

These results provide justification for the 
stability of relativistic one boson exchange 
models that usually exclude matter loops but 
may include Z-diagrams of all orders. Our 
argument cannot be easily extended to sym- 
metric (p 3 theories where it is impossible to 
make a clear distinction between Z-diagrams 
and loops. 

VI. CONCLUSIONS 

In this paper we have given a summary 
results for scalar interactions obtained with 
the use of the FSR representation . The FSR 
approach uses a covariant path integral rep- 
resentation for the trajectories of particles. 
Reduction of field theoretical path integrals 
to path integrals involving particle trajecto- 
ries reduces the dimensionality of the prob- 
lem and the associated computational cost. 

Applications of the FSR approach to 1 and 
2-body problems in particular shows that un- 
controlled approximations in field theory may 
lead to significant deviations from the correct 
result. Our results indicate that use of the 
Bethe-Salpeter equation in ladder approxi- 
mation to solve the 2-body bound state prob- 
lem is a poor approximation. For the scalar 
theories examined here, a better approxima- 
tion to the 2-body problem is obtained us- 
ing the Gross equation in ladder approxima- 
tion. Similarly, use of the rainbow approxi- 
mation for the 1-body problem gives a poorer 
result that simply calculating the self energy 
to second order! In all of these cases, the ex- 
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planation for these results seems to be that 
the crossed diagrams (such as crossed lad- 
ders) play an essential role, canceling contri- 
butions from higher order ladder or rainbow 
diagrams. 
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